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Abstract 

We consider the quasiconformal dilatation of projective transformations of 
the real projective plane. For non-affine transformations, the contour lines of 
dilatation form a hyperbolic pencil of circles, and these are the only circles that 
are mapped to circles. We apply this result to analyze the dilatation of the cir- 
cumcircle preserving piecewise projective interpolation between discretely con¬ 
formally equivalent triangulations. We show that another interpolation scheme, 
angle bisector preserving piecewise projective interpolation, is in a sense opti¬ 
mal with respect to dilatation. These two interpolation schemes belong to a 
one-parameter family. 

30C62, 52C26 

1 Overview and motivation 

The deviation of a function from being conformal is measured by its quasiconformal 
dilatation. (Standard ways to quantify the local deviation from conformality will 
be reviewed in Section 2.) This paper is about the quasiconformal dilatation of 
real projective transformations of the plane. In Section 3 we show that the contour 
lines of dilatation form a hyperbolic pencil of circles (Theorem 3.1). These are the 
only circles that are mapped to circles (Theorem 3.4). Although the statements 
are elementary and the proofs are straightforward, these results seem to be new 
It may well be that the dilatation of a real projective map was never considered 
before. We were motivated by the theory of discrete conformal maps that is based 
on the following definition of discrete conformal equivalence of triangle meshes. 
(Quasiconformal dilatation also plays an important role in at least one other theory 
of discrete conformality, circle packing [8].) 

1.1 Definition and Theorem ([3, 5]). Two triangulated surfaces are considered 
discretely conformally equivalent, if the triangulations are combinatorially equivalent 
and if one (and hence all) of the following equivalent conditions hold (see Figure 1 
for notation). 

(i) The edge lengths and I^ of corresponding edges are related by 



for some logarithmic scale factors iq associated to the vertices, 
(ii) For interior edges ij, the length cross ratios are equal, that is, 

^mj^ki imj^ki 
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Figure 1: Notation for the definition of discrete conformal equivalence. 


(iii) The circumcircle preserving projective maps that map a triangle of one triangu¬ 
lation to the corresponding triangle of the other triangulation, and the respec¬ 
tive circumcircles onto each other, fit together continuously across edges. 

The original definition (i) is due to Feng Luo [5]. For the equivalent charac¬ 
terizations (ii) and (iii), see [3]. Characterization (iii) means that discretely con¬ 
formally equivalent triangle meshes allow not only the usual piecewise linear (PL) 
interpolation (which works for any two combinatorially equivalent triangulations) 
but also circumcircle preserving piecewise projective (CPP) interpolation (which is not 
continuous across edges unless two triangulations are discretely conformally equiv¬ 
alent) . This motivated the following definition. 

1.2 Definition ([3]). A simplicial continuous map between triangulated surfaces is 
a discrete conformal map if the restriction to any triangle is a circumcircle preserving 
projective map onto the image triangle. 

Figure 2 shows visualizations of PL and CPP interpolations. The CPP interpola¬ 
tions clearly look better. (This is an important point for applications in computer 
graphics [7].) What is the reason? This question led us to investigate the quasicon- 
formal dilatation of projective transformations. We wanted to check the hypothesis 
that the CPP interpolations had lower dilatation. Behind this hypothesis was the 
non-mathematical hypothesis that lower dilatation was the reason why CPP inter¬ 
polation looks better; however, see Remark 6.4 in Section 6. 

As we show in Section 4, the dilatation of CPP interpolation is indeed pointwise 
less than or equal to the dilatation of PL interpolation, but the maximal dilatations 
per triangle are equal (Corollary 4.2). 

In Section 5, we show that another interpolation scheme, mapping angle bisec¬ 
tors to angle bisectors, is in a sense optimal with respect to dilatation (Theorem 5.1). 
Like CPP interpolation, angle bisector preserving piecewise projective (APP) inter¬ 
polation is continuous across edges if and only if the triangulations are discretely 
conformally equivalent (Theorem 5.3). 

The interpolations schemes PL, CPI? and APP are in fact members of a continuous 
one-parameter family (see Section 6). We do not know any interesting geometric 
characterization for any other member of this family. 


2 Dilatation and eccentricity 

In this section, we review the definitions of dilatation and eccentricity of real dif¬ 
ferentiable maps. The definitions are standard (the provided references are only 
meant as examples), but our perspective is slightly unusual because (in Section 3) 
we are interested in maps that become singular and orientation reversing. (Real 
projective transformations are generally not quasiconformal functions on C.) 
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Figure 2: Piecewise linear (PL) vs circumcircle preserving piecewise projective 
(CPP) interpolation. Top row [7]: A triangulated surface in R 3 resembling a cat 
head is mapped to a conformally equivalent planar triangulation (not shown) by 
PL and CPP interpolation. The maps are visualized by pulling back a checkerboard 
pattern in the plane to the cat head. Some notable differences are highlighted. 
Middle row [3]: A triangulated planar region is mapped to a conformally equiva¬ 
lent triangulation of a rectangle (bottom) by PL and CPP interpolation. 
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In the following, let U c C be open and let f : U —> C be a real differentiable 
map. (We identify C and M 2 as euclidean vector spaces.) The function / maps a 
small circle of radius s around z eU approximately to an ellipse with major semi¬ 
axis X 1 e and minor semi-axis X 2 e, where 


X 1 > A 2 > 0 


are the singular values of the real derivative df z of f at 2 . 

The quotient of singular values is usually called the dilatation [1, 2], or, to be 
more specific, the quasiconformal dilatation [6]. However, since we are interested 
in maps that may become singular and orientation reversing, we define the [signed] 
dilatation Df(z) £ R U { 00 } of f at the point z eU by 

Ai 

D/(z) = ±—, 

where the sign is chosen according to whether f is orientation preserving (+) or 
reversing (— ) at 2 . If the derivative df is singular at z, then X 2 = 0 and we define 
Df(z) = 00 . But we will assume that df vanishes nowhere, so the dilatation is well 
defined everywhere. Note the following properties of the dilatation: 

• \D f \>l, 

• Df =± 1 where / is conformal or anticonformal, respectively, 

• D / (z) = D / -i(/(z)). 

The Beltrami coefficient 

M/0) = / i // 2 (1) 

is defined in terms of the Wirtinger derivatives 

fz = ffx - ify), fz = ffx + ify )• 

The modulus of the Beltrami coefficient is called the eccentricity of f [1]. The ec¬ 
centricity \(Jif | is related to the dilatation Df by 

Df-l 
D f + 1' 

Note that \(j,f\ e [0, 00 ], and 

JO where f z = 0, that is, where f is conformal, 

* ) 00 where f z = 0, that is, where f is anticonformal, 

\[j,f | = 1 where det df = 0, 


and 

IjU/OOl = lM/-i(/(z))|. (2) 

2.1 Remark (on terminology), (i) A function is called quasiconformal if it is ori¬ 
entation preserving and its dilatation is bounded. In phrases like “quasiconformal 
dilatation”, the adjective “quasiconformal” is used somewhat sloppily to mean “be¬ 
longing to the theory of quasiconformal functions”, (ii) The above definition of 
“eccentricity” has fallen into disuse. The Beltrami coefficient [if is also called com¬ 
plex dilatation [2], although complex eccentricity would make more sense. 


4 




Figure 3: The contour lines of |pif | are a hyperbolic pencil of circles (Theorem 3.1). 


3 Dilatation of a projective map, and circles mapped 
to circles 

We are interested in the dilatation of a projective map 

f : RP 2 RP 2 , [x] [Ax], A e GL 3 (R), 

where we identify the complex plane C with the real projective plane RP 2 without 
the line t^ = {[x] | x 3 = 0} via the map 

z = x + iy «— >[y]. 

This provides the conformal structure on RP 2 \ t ^. We consider the line t ^ as 
the line at infinity. The map f is [real] affine if it maps £ 00 to The map f 
is a similarity transformation if it is real affine and complex affine on C, that is, 
% »-> az + b with a £ C \ {0}, b e C. 

The eccentricity of affine maps is constant, and identically 0 for similarity trans¬ 
formations. The following theorem treats the interesting case (see Figure 3). 

3.1 Theorem. If the projective transformation f : RP 2 —> RP 2 is not affine then: 

(i) The contour lines of \pif \ are a hyperbolic pencil of circles. 

(ii) The function \pif\ is an affine parameter for this 1-parameter family of circles 
that assigns the values 0 and oo to the circles that degenerate to points and the 
value 1 to the preimage of the line at infinity under f. 

(iii) The function f maps this hyperbolic pencil of circles to another hyperbolic pencil 
of circles 

The following proof relies on direct calculations. It would be nice to have a 
more conceptual argument. 

Proof In terms of the affine parameter % on RP 2 \ f oo, the projective map / may 
be written as 

az + bz + c 

z ■—>-rz-, a,b,c,peC, q £ R, (3) 

pz + pz + q 
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where p ^ 0 because f is not affine. In terms of the coefficients a, b,c,p,q , the 
matrix representing f is 


A = 


^Q-i + bi 

a 2 + b 2 

V 2 pi 


—a 2 + b 2 
a i ~ b 1 
-2 p 2 



where indices 1 and 2 to denote real and imaginary parts, respectively. 

From the definition (1) of p, using the representation (3) for /, one obtains 
immediately 


M/0) = 


az + /3 
—az + y ’ 


where 


{ a = bp — ap 
I5 = bq-cp , 
T = aq — cp 


SO 

az + f3 
—az + j 

As a tedious but straightforward calculation shows. 


IM/OOI = 


det Q=pdetA#0 3 

so \pf(z)\ = \M(z)\, where M is a Mobius transformation. Parts (i) and (ii) of the 
theorem follow easily. 

To see part (iii), note that the inverse map / -1 is again a projective map which is 
not affine. Therefore the contour lines of \pf-i\ are also the circles of a hyperbolic 
pencil of circles. Because f and / _1 have the same dilatation at corresponding 
points (equation (2)), f maps circles with constant \pf\ to circles with constant 
\Pf-i\. This proves (iii). □ 

3.2 Corollary. If f is orientation preserving on a triangle T then the function \pf\ 
attains the maximum 

max||U/(z)| 

at a vertex. 


Indeed, in the open half-plane where f is orientation preserving, the sublevel 
sets of \pf | are strictly convex (being disks), and every point is on the boundary 
of its sublevel set. Hence, the maximum cannot be attained in the interior of the 
triangle or in the relative interior of a side. 

3.3 Remark. Here and throughout this article, “triangle” shall mean “closed trian¬ 
gular region”. 

Theorem 3.1 suggests the following question: Which circles are mapped to cir¬ 
cles by a projective map /? Obviously, a similarity transformation maps all circle to 
circles, and an affine transformation that is not a similarity maps all circle to ellipses 
that are not circles. The following theorem treats the remaining case. 

3.4 Theorem. If the projective transformation f : RP 2 —> RP 2 is not affine then it 
maps precisely the circles of a hyperbolic pencil to circles. 
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Proof. This proof relies on the classical characterization of circles in terms of their 
complex intersection points with the line at infinity. Consider the real projective 
plane RP 2 as the set of points with real homogeneous coordinates in the complex 
projective plane CP 2 . In CP 2 , every conic section intersects the line at infinity twice 
(counting multiplicity). Circles are the conics intersecting the line at infinity in the 
imaginary circle points 



The circles that f maps to circles are the conics containing K, K, / _1 (fC), and 
/ _1 (K). Because f is real and not affine, these four points are in general position: 
No three of them are contained in a line, which would have to be the line at infinity 
or its preimage under f . But these lines intersect in a real point. Hence the conics 
containing the four points form a pencil of conics. Since the four points are two 
pairs of complex conjugates, the pencil is a hyperbolic pencil of real conics. □ 


4 The circumcircle preserving projective map 

The following theorem characterizes the maximal value of \pif\ for a circumcircle 
preserving projective map on a triangle. 

4.1 Theorem. If the projective map f maps a triangle T with vertices A, B, C onto 
a triangle f with vertices A, B, C, and the circumcircle of T to the circumcircle of f, 
then 

\Hf(A)\ = Im/(B)I = lM/(C)l = IjUfcl, (4) 

where h is the real affine map from T onto t (whose eccentricity \p h \ is constant). 

Together with the results of the previous section, Theorem 4.1 implies the fol¬ 
lowing corollary. 

4.2 Corollary. If a projective map f satisfying the assumptions of Theorem 4.1 is 
orientation preserving on T, then 


lM/(P)l < lMh(-P)l; 

for all points PET, and 

max|M/(P)| = \n h \. 

4.3 Remark. The projective map f of Theorem 4.1 is either orientation preserving 
or orientation reversing on all of T. Since f maps circumcircle to circumcircle, 
neither nor / _ 1 (^oo) intersects T. 

Proof of Theorem 4.1. Since the circumcircle is mapped to the circumcircle, it is a 
contour line of \pif | (Theorems 3.1 and 3.4). This proves the first two equalities of 
equation (4). It remains to show the last equality involving the eccentricity of the 
affine map. We will avoid any calculation of derivatives. 

Given three lines intersecting in one point P and three lines 

intersecting in one point P, there exists an affine map F mapping £j to (j = 
1,2,3), and this map is unique up to post-composition with a homothety centered 
at P (or, which amounts to the same, pre-composition with a homothety centered 
at P). The value of p F is therefore uniquely determined by the lines. 
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Figure 4: The affine map h from T onto T. 




Figure 5: The circumcircle preserving map f from T onto f. 

The affine map h maps the point A to the point A and the lines £ 1 = AB,l 2 = AC, 
and the line i 3 parallel to CB through A to the lines 1 1 = AB, l 2 — AC, and the line 
l 3 parallel to CB through A (see Figure 4). 

The circumcircle preserving projective map f also maps i 1 to £ 1 and i 2 to l 2 , and 
it maps the tangent 1' 3 to the circumcircle at A to the tangent 1' 3 to the circumcircle 
at A (see Figure 5). The affine approximation of f at A (that is, in affine coordinates, 
the first order Taylor polynomial) also maps these lines to the same lines. Consider¬ 
ing the angles at A and A in Figures 4 and 5, one finds that the affine approximation 
of / at A is equal to h up to composition with similarity transformations. Thus, 
|jU/(A)| = K|. □ 

5 The angle bisector preserving projective map 

Which projective transformation between two given triangles minimizes the maxi¬ 
mal dilatation? Since the maximal dilatation is attained at a vertex (Corollary 3.2), 
it is enough to minimize the maximal dilatation at the vertices. As it turns out, there 
is a unique projective transformation that simultaneously minimizes the dilatation 
at all three vertices: 

5.1 Theorem. Of all projective maps that map a triangle T cC with vertices A, B, C 
onto a triangle f c C with vertices A, B, C, which are labeled in the same orientation 
so that the maps are orientation preserving, the one that maps the angle bisectors to 
the angle bisectors simultaneously minimizes |jU(A)|, |ju(£>)|, and |ju(C)|. That is, the 
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Figure 6: The angle bisector preserving projective transformation maps incircle cen¬ 
ter to incircle center, but it does in general not map incircle to incircle. 



Figure 7: Angle bisector theorem: | = |. 


angle bisector preserving projective transformation f satisfies 

lM/(A)| < |ju g (A)|, ||U/(B)I < \n g (B)\, \n f (C)\ < |jU g (C)| 

for all projective transformations g with g(A) = A, g(B) = B, g(C) = C, g(T) = f. 

This theorem follows immediately from the following Lemma, whose proof we 
leave to the reader. (All arguments we know rely ultimately on one or another more 
or less direct calculation.) 

5.2 Lemma. Of all linear maps in SL 2 (R) that map two one-dimensional subspaces 
Rv, R w onto two one-dimensional subspaces Rv, Rw, the map f e SL 2 (R) has the least 
dilatation if and only if it maps one , and hence both , of the angle bisectors R( ^ ± ^ ) 

to the corresponding angle bisector r(j^ ± j^j). (This determines f uniquely up to 
sign.) 

The characterization of discrete conformal equivalence in terms of projective 
maps (Definition and Theorem 1.1 (iii)) remains true if “circumcircle preserving” is 
replaced by “angle bisector preserving”: 

5.3 Theorem. The angle bisector preserving projective maps between corresponding 
triangles of combinatorially equivalent triangulations fit together continuously across 
edges if and only if the triangulations are discretely conformally equivalent. 

Proof. This follows immediately from the elementary angle bisector theorem (Fig¬ 
ure 7) and the characterization of discrete conformal equivalence in terms of length 
cross ratios (Definition and Theorem 1.1 (ii)). □ 

6 A one-parameter family of piecewise projective in¬ 
terpolations 

The angle bisector preserving projective transformation maps incircle center to in¬ 
circle center. The incircle center in a triangle ABC has barycentric coordinates 
[a, b, c], where a, b, c denote the lengths of opposite sides. 
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Figure 8: The circumcircle preserving projective transformation mapping ABC to 
A'B'C' maps the point S to the point S'. The points S and S' are called the symme- 
dian, Lemoine, or Grebe point of the triangles ABC and A'B'C', respectively. (They 
are the Gergonne points of the triangles formed by the tangents.) 


The circumcircle preserving projective transformation maps symmedian point to 
symmedian point (see Figure 8). The symmedian point has barycentric coordinates 
[a 2 , b 2 ,c 2 ]. 

In the affine case, the barycenter is mapped to the barycenter. Its barycentric 
coordinates are [1,1,1]. 

6.1 Definition. For t e R, let the exponent-1-center in a triangle ABC be the point 
with barycentric coordinates [a 1 , b l ,c l ], where a, b, c are the lengths of opposite 
sides. 

6.2 Remark. For parameter values different from t = 0 (barycenter), t = 1 (incircle 
center), and t = 2 (symmedian point), the exponent- 1 -centers of a triangle seem 
to be fairly esoteric triangle centers. For example, the values t = 3,4 correspond 
to triangle centers X(31) and X(32) in Kimberling’s list [4], and the values t = 
-1,-2 correspond to X(75) and X{7 6). If there are any other exponent- 1 -centers 
in Kimberling’s list, they have indices greater than 300. 

The proof of Theorem 5.3 obviously generalizes to “exponent- 1 -center preserv¬ 
ing projective maps”: 

6.3 Theorem. For any teM, the projective maps between corresponding triangles of 
combinatorially equivalent triangulations that map exponent-t-centers to exponent-1- 
centersfit together continuously across edges if and for t 0 only if the triangulations 
are discretely conformally equivalent. 

This leads to a one-parameter family of exponent- 1 -center preserving piecewise 
projective interpolation schemes for discretely conformally equivalent triangula¬ 
tions. The parameter value t = 0 corresponds to piecewise linear interpolation, 
t = 1 to angle bisector preserving, and t = 2 to circumcircle preserving piecewise 
projective interpolation. Figure 9 visualizes these interpolation schemes for some 
values of t using the same example as Figure 2, middle. 

6.4 Remark. In our eyes, circumcircle preserving piecewise projective interpolation 
(t = 2, CPP) looks slightly better than angle bisector preserving interpolation (t = 
1, APP). We have found this in other examples as well. Since APP has in general 
lower maximal dilatation per triangle, this indicates that low dilatation is not what 
makes CPP interpolation look better. We do not know which mathematical property 
of CPP interpolation accounts for the nicer appearance. 
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Figure 9: One-parameter family of piecewise projective interpolations. 
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